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Abstract 



The path integral approach is used for the calculation of the correlation functions of the 
XY Heisenberg chain. The obtained answers for the two-point correlators of the XX 
magnet are of the determinantal form and are interpreted in terms of the generating 
functions for the random turns vicious walkers. 
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1 Introduction 



The problem of enumeration of paths of vicious walkers on one-dimensional lattice was 
formulated by M. Fisher [1] and since then continues to attract much attention (see refs. 
in [2,3]). The walkers are called 'vicious' because they annihilate each other at the same 
lattice site, and their trajectories are thus non-intersecting. Similar problems appear in the 
theory of domain walls [4] , directed percolation [5] , self-organized criticality [6] , and polymer 
theory [7] . It has been proposed Ref . [2] to use the XX Heisenberg chain to enumerate the 
paths of the random turns vicious walkers. 

Approach based on path integration was developed in Refs. [8, 9] to calculate thermal 
correlation functions of the XY Heisenberg magnet. Dependence of the integration variables 
on the imaginary time is defined by special quasi-periodicity conditions. In the present paper, 
this method is used for the calculation of the two-point correlation functions of the XX model 
and the interpretation of the obtained answer in terms of the generating functions of the 
random turns vicious walkers is given. 



2 The problem 

The Hamiltonian of the periodic XY Heisenberg chain of "length" M (M is chosen to be 
even) in transverse magnetic field /i > is: 

M 



n,m=l 

^ M ^ M 
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Here is z-component of the total spin operator, and the entries of the so-called hopping 
matrix Anm {s — ±) are: 

2^nm = ^\n-m\,l + S(5|„_m|,M-l , (2) 

where 6n,i is the Kronecker symbol. The Pauli matrices = (1/2) (a^ ± ia^) and a^, 
where n E A4 = {!,..., M}, satisfy the commutation relations: [cr^, a^] = Skiaf and 
a^^ — ±2SkiO'^- The periodic boundary condition reads: ct^+a^ = cr^, Vn. The 
Hamiltonian H (1), taken at 7 = (the case of XX magnet), commutes with S^. 
Time-dependent thermal correlation functions are defined as follows: 

G„,(m, t) ^ Z-'Tr «„(0) af (t) e'^^) , Z^Tr {e'^^) , (3) 

where ^^(t) = e**^ cif e^**-'^, f3 = 1/T is inverse temperature, and t is time. This correlator 
may be rewritten in terms of the canonical lattice Fermi fields q, cj, where i, j E Ai, hj 
means of the Jordan- Wigner map: 

n—l n—l 
= (n ^l) ^'^ ' ^" = (n ^j) ' neM, 

where — 1 — ^CjCj. The periodic conditions for the spin operators result in the boundary 
conditions for the fermions: 

cm+1 = (-l)'^ci, cj^^+i = cl(-l)-^, (4) 
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where N = X^^li cJjC„ is the operator of the total number of particles. In the fermionic 
representation, H (1) will take a form H = H+P+ + R-p-, where = (V2)(I ± (-l)"^) 
are projectors [8]. The operators are of identical form with s — ± pointing out a 
correspondence between these operators and appropriate specification of the conditions (4) : 

cm+1 = "-sci, c^^^i^i = —sc\. 

Equation (3) for the 2;-components of spins, for instance, becomes: 

G,,im,t) = 1 - 2Z-ilV(c;+„Q+^e-^^) - 2 ^-^IV (cf q e"^^) 

+ 4 Z-'Tt {cl^^ Q+„, e^*^ 4 Q e-(/5+'*)^) . (5) 
To evaluate (5), it is convenient to consider the generating functional: 

g = g{S, T\ II, v) = Z-^IV (e^ e-^"" e"^^) , (6) 

where yU, v are the complex parameters, /i + z/ = /5. Two operators, S = Sc and T = c^Tc, 
are defined through the matrices S = diag [Si, S2, ■ ■ ■ , Sm}, T = diag |Ti, T2, . . . , Tm}- 
For instance, the last term in R.H.S. of (5) is obtained from (6) in the following way: 



T„,Vn 
II, V — > —it, (3 + it 



As a result, we express the trace in R.H.S. of (6) in the form [8]: 

IV (e^ e-'^^ e""^^) = ^ (g^ Z+ + + Z+ - g^ Z^) , (7) 

where 

= Tr(e^e-''^%^e-'^^*), 
= Tr(e^e-^^%^(-l)^e-'^^*), 
and Z± = Tr(e-'3^*), = Tr ((-l)^e-^^*). 



3 The path integral 

We use the coherent states 1^;) = exp(c^2;)|0) and (^z*\ = (0| exp(2;*c) generated from the 
Fock vacuum |0), Cfc|0) = 0, V/c. We use the short-hand notations for the M-component 
objects, say, z* = {zl, . . . , z\j) and z = {zi, . . . , zm) formed by the independent Grassmann 
parameters Zk,zl. [k G M.). Besides, Xlfli '^I-^fc = ^^^^ Ilfcli = Then, we shall 

represent [9] the trace of the operator in g^ (8) by means of the Grassmann integration 
over d^;, dz*: 

g^Z^ = J dzdz*e'''{z*\e^e-''^^ e^e-^"^\z). (9) 

For the sake of simplicity we shall consider the XX model only and take those that 
correspond to iJ (1) at 7 = 0. 

To represent R.H.S. of (9) as the path integral, we first introduce new coherent states 
\x{I)), (a;*(/)|, where 2Lx M independent Grassmann parameters are arranged in the form 
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of 2L "vectors" x*{I), x{I) (/ e {1, . . . , L}). It allows to insert L times the decompositions 
of unity 

dx%I)dx{I) exp{-x*{I)x{I)) \x{I)){x*{I)\ 



into R.H.S. of (9). We define then the additional variables satisfying the quasi-periodicity 
conditions: 

-Ex{0) ^ x{L + l) = z, -x*{L + l) ^ x*{0)E-^ = z* . (10) 

Here, E = e"^^"*" with the matrices if^ expressed [9] through the hopping matrices 
(2): = - AW + hi, where / is a unit M x M matrix. The described procedure allows to 
pass in the limit L — oo from {L + l)-fold integration to the continuous one over "infinite" 
product of the measures dx*{T)dx{T) on a space of trajectories x*{t), x{t), where r e M: 



g^Z^^ J e^dX*dXY[dx*{T)dx{T). 



The integration over the auxiliary Grassmann variables A*, A guarantees the fulfilment of the 
continuous version of the constraints (10). The action functional is 5* = J L{T)dT, where 
L{t) is the Lagrangian: 

L(t) = x*{t) - x{t) + r{T)x{T) + x*{t)J{t) , 

J*(r) = \*{5{t)I + 5{r-u)E-^) , J{t) = {S{t)T + 5{r-u)E)\. 

The 5-functions reduce r to the segment [0, (3]. The stationary phase requirements 6S/6x* — 
0, 6S/6x = yield the regularized answer [9]: 

- det / + ^ ^- , 

The remainder correlators Ga(,(m,t) (3) (with a, 6 e {+, — }) are obtained analogously. 



4 Random walks 

The evolution of the states obtained by selective hipping of the spins governed by the XX 
Hamiltonian Hq (1) is related to a model of a random turns vicious walkers [2,3]. Indeed, 
let us consider the following average over the ferromagnetic state vectors (ft" | , | ft") : 

F,,KA) = (^k;e-^^Vrlt), (11) 

where | '(|) = ])ni i-©-) all spins arc up, and A is an "evolution" parameter. Spin up 

(or down) corresponds to empty (or occupied) site. Differentiating Fj.i{X) (11) and applying 
the commutator [Hq, ct^], we obtain the differential-difference equation {master equation): 

^ Fr,iW = I {Fj+i;iW + Fj-i;iW) ■ (12) 

The average Fj-i may be considered as the generating function of paths made by a random 
walker travelling from l^^ to j*^ site. Really, its K-ih derivative has the form 

A=0 ni,...,nK-i 



d^ 
d\^ 



4 



A single step to one of the nearest sites is prescribed by the hopping matrix (2) with s = +. 
After K steps, each path connecting Z*^ and j*^ sites contributes into the sum. The A'"-point 
correlation function [N < M), 



(13) 



enumerates the nests of the lattice paths of N random turns vicious walkers being initially 
located at the positions li > I2 > ■ ■ ■ > In and, eventually, at ji > 32 > ■ ■ ■ > 3n- It is 
expressed in the form [2] : 



l<r,s<N ■ 



(14) 



The ground and the excited states of the XX chain at h — with the total spin equal 
to (M/2) — N are decomposed over a basis of states (^fi'^i^ ■ ■ ■ with spins flipped 

[10]. Therefore, the trace -Fm+i:i(A) = Tr (o"^_,_^e^'*'-^"o'f) is a linear combination of the 
generating functions (13) describing the evolution of A^ + 1 random turns walkers. The 
initial and the final positions of one of them are fixed at /i = 1 and ji = m + 1, respectively, 
while for the rest ones these positions are random. In the thermodynamic limit, the number 
of the virtual walkers tends to infinity. We apply the procedure described in 3 to calculation 
of -F'm+i;i(A) in the limit when M and A^ are large enough. In this limit, the contribution 
of the terms with the subindex 'B' become, with regard at (8), negligible in (7). We thus 
obtain: 



F^+i;i(A) = tr (e-^-^ ei,^+i) 



M 



) 



a=0 



det(/+W^) 



.(e--.,„,.)-l.(^) 



where ei^m+i = {Si^nSm.+i,i)i<n,i<M, and the matrix is used instead of since s can 
be taken zero at large enough M. The traces of A-dependent M x M inatrices I4m and Vm 
are given below. Differential equation analogous to (12) is fulfilled by F^+i;i(A). At large 
separation m it takes the form: 



^F^+i; i(A) = i (F^; i(A) + i(A)) - Tr {Ho a+^. 



(15) 



We expand formally F^+i; i(A) with respect to Um and obtain the answer in two lowest orders 
as follows: 



^m+l; l(A) 



m+l;l(A) XtlUm - J^^'^^rn, 



Fm+1; 1 (A) + F, 
(M - 2m)Fi;i(A), 

m 

Fm+1; 1 (2A) - 2 5^F„+i;;(A)F;;i(A). 



(16) 



1=1 



Although M and m are chosen to be large in this expansion, the ratio m/M is kept bounded. 
In each order the master equation (15) is fulfilled by (16). The contribution of the second 
order can be re-expressed through the two-point functions Fm+i,/;i,i(A) (see (13), (14)). Thus, 
summation over intermediate positions (of a virtual walker located at site) arises in the 
second order. A similar picture is expected in next orders. 
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